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Abstract. We review the interaction in intermediate and high mass stars between their evolution 
and magnetic and chemical properties. We describe the theory of Ap-star 'fossil' fields, before 
touching on the expected secular diffusive processes which give rise to evolution of the field. We 
then present recent results from a spectropolarimetric survey of Herbig Ae/Be stars, showing that 
magnetic fields of the kind seen on the main-sequence already exist during the pre-main sequence 
phase, in agreement with fossil field theory, and that the origin of the slow rotation of Ap/Bp 
stars also lies early in the pre-main sequence evolution; we also present results confirming a lack 
of stars with fields below a few hundred gauss. We then seek which macroscopic motions compete 
with atomic diffusion in determining the surface abundances of AmFm stars. While turbulent 
transport and mass loss, in competition with atomic diffusion, are both able to explain observed 
surface abundances, the interior abundance distribution is different enough to potentially lead 
to a test using asterosismology. Finally we review progress on the turbulence-driving and mixing 
processes in stellar radiative zones. 

Keywords, stars: magnetic fields - stars: chemically peculiar - stars: evolution - (magnetohy- 
drodynamics) MHD - stars: rotation - stars: intermediate mass - stars: pre-main sequence 



1. The interesting lives of middle main sequence stars 

Observations reveal that around 5-10% of intermediate-mass main-sequence (MS) stars 
(1.6 - 8M Q ) have magnetic fields of strength 300 G - 30 kG (see Donati & Landstreet 
2009 for a recent review of the so-called 'chemically peculiar' (CP), or Ap and Bp stars). 
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The fields are large-scale and apparently static, i.e. they do not evolve at all over periods 
of decades, and their host stars rotate much more slowly than normal A stars. 

It is worthwhile to look at magnetic middle MS stars in the broader context of stellar 
evolution. Stars form by gravitational collapse on a timescale of a few 0.1 Myr, set by 
the free-fall time. At a size of a few 100 AU, angular momentum leads to slower growth 
through quasi-static disk accretion. The proto-star shrinks towards the MS as a T Tau 
star, and later, if massive enough, as a Herbig AeBc star. During early stages of collapse 
much magnetic flux may be retained, or lost via ambipolar diffusion, but it is not clear 
how to retain flux whilst the star is almost fully convective, or how the field of some 
Herbig AeBe stars might be related to the earlier T Tau dynamo activity phase. 

On the MS, surface chemistry changes occur as a result of gravitational and radiatively 
driven atomic diffusion, which must compete with convection, meridional circulation and 
turbulent diffusion, mass loss from the surface, the magnetic field (especially in the upper 
atmosphere and beyond) , and possibly even accretion. The complications of this multi- 
dimensional competition probably explain why some stars with fossil fields are marked 
by very distinctive chemical peculiarities (hence 'CP' stars) while other, hotter magnetic 
stars are not, but the details of the effects are still far from clear. However, the way 
in which different effects predominate with changes in stellar parameters such as mass, 
age and rotation rate make these stars excellent laboratories for the study of the various 
important internal physical processes. 

Magnetic fields have been known in peculiar B stars up to about 8 M Q for decades, 
but are now being found in a (still small) fraction of stars well above this mass, mostly 
without chemical peculiarities (e.g. Bouret et al. 2008 and refs. therein). These fields have 
significant effects on the strong stellar winds, and in extreme cases wind gas is trapped 
in closed field loops (Landstreet & Borra 1978; Babel & Montmerle 1997); this effect has 
been beautifully modelled in detail by Townsend et al. (2007). 

The host stars of these fossil fields straddle the boundary between low mass stars which 
have greatly increased luminosity as giants, and high-mass stars that evolve at almost 
constant luminosity as supergiants. Note, that during the MS phase, R increases by a 
factor of 2, and T e g decreases by 30%. The magnetic fields also pose major puzzles. Why 
are no fields, dynamo or fossil, found in MS F3 - F5 stars? Why do fields in MS Ap 
stars appear to decay on a time scale short compared to the Cowling time but still a 
significant fraction of the MS lifetime (Landstreet et al. 2008)? How and when do strong 
magnetic fields influence rotation rates? What happens to fields as the host stars evolve 
to the giant phase - can the fossil field anchor to stable layers deep in the star and persist 
in spite of the deep convection? Can fossil fields survive the giant phases to the end, to 
become the fossil fields of white dwarfs, or do these have a different origin, perhaps in a 
close binary system? 

2. Fossil fields in Ap stars and degenerate stars 

Over the 60 years since magnetic fields were discovered in Ap stars the subject has ac- 
companied the development of magnetohydrodynamics, as theoretical explanations have 
been sought to explain the properties of the observed 'oblique rotators' - large-scale, 
static fields. Two theories have been proposed. In core-dynamo theory the field is gener- 
ated in the core by differential rotation and convection, as is known to work in the solar 
envelope, and rises through the radiative zone to the surface - and herein lies a severe 
drawback (another being the lack of positive correlation between field strength and ro- 
tation speed): There are diffusive mechanisms to cause this rise, but the timescales are 
apparently too long. According to fossil theory, what we see is a remnant from either 
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the original molecular cloud or a protostellar dynamo - when the star settles down, the 
field relaxes into a stable equilibrium, which then evolves only via the same very slow 
diffusive processes over timescales longer than the MS lifetime. 

Theoretical work on fossil fields has focussed on demonstrating that such stable equi- 
libria are in fact possible. This is problematic, as it can be shown that the simplest 
magnetic field configurations are always unstable (Goossens et al. 1981). Generally stud- 
ies have looked only at axisymmetric fields: Tayler (1973) showed that all purely toroidal 
fields are prone to interchange (axisymmetric) and kink (non-axisymmetric) instabilities, 
and Markey & Tayler (1973) and Wright (1973) showed that purely poloidal fields are 
also unstable. Therefore it was concluded that a stable axisymmetric field (if it exists) 
must indeed have both poloidal and toroidal components. Now, such a field can be ex- 
pressed as the sum of a poloidal and a toroidal component, each completely described 
by a single scalar function. The Lorentz force for an axisymmetric field cannot have an 
azimuthal (cf>) component, since no counterpart exists in the hydrostatic fluid forces that 
can balance it. From this it can be shown that the two functions defining the poloidal 
and toroidal parts must be related. One consequence of this is that the toroidal field 
is entirely confined within the poloidal field lines that close inside the star. Simulations 
have now demonstrated that in stably stratified stars, random initial magnetic fields can 
evolve into roughly axisymmetric configurations with both poloidal and toroidal compo- 
nents, which then remain stable over a diffusive timescale (Braithwaite & Spruit 2004). 
It would be desirable to provide an analytic proof for the existence of such fields, and to 
understand how the poloidal and toroidal components help to stabilise each other. Using 
the energy principle developed by Bernstein et al. (1958), we hope to be able to show that 
instabilities present for toroidal fields can be eliminated by adding a poloidal component, 
and vice versa (Akgiin & Reisenegger, in prep.) This could also allow us to determine the 
range of the ratio of the poloidal and toroidal components. Some progress towards that 
end (using a mixture of analytic and numerical methods) has been made by Braithwaite 
(2009), where it is found that the toroidal component can be either much stronger than, 
or of comparable strength to, but not much weaker than, the poloidal component. Mean- 
while, non-axisymmetric equilibria have been found in simulations (Braithwaite 2008). 
Even in the most strongly magnetised stars, the Lorentz force is still typically a million 
times weaker than hydrostatic forces due to pressure and gravity. Therefore, a small 
perturbation in the background equilibrium is sufficient to balance the Lorentz force. In 
radiative envelopes of massive stars and in degenerate stellar interiors, matter is non- 
barotropic and stably stratified, which hinders radial displacements - this allows a wider 
range of magnetic field structures than found in barotropic fluids (Reisenegger 2009). In 
Braithwaite et al. (in prep.) we show that a star composed of matter with a barotropic 
equation of state is unlikely to retain a fossil field even if equilibria are available, since 
magnetised regions rise buoyantly and are lost through the surface on a timescale com- 
parable to the reconnection timescale. We also show that analogous equilibria can exist 
outside of stars, for instance inside radio bubbles inflated by AGN - see Fig. [TJ 

It is informative to look at the process of relaxation of an arbitrary field into a stable 
equilibrium. Equilibria have special topological properities, such as the zero-gravity case 
where equilibrium requires that VP = (l/47r)(VxB) xB and, since the Lorentz force is 
perpendicular to the field, the field lines reside in isobaric magnetic surfaces. Relaxation 
therefore requires topological change, i.e. magnetic reconnection. Magnetic helicity H = 
J B • A dV where B = V x A is approximately conserved during relaxation, enabling us 
to predict the energy (and field strength) of the final equilibrium if we know the helicity 
of the initial turbulent field, since the energy and helicity of the equilibrium are related 
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Figure 1. Field lines of 
two magnetohydrostatic equi- 
libria. The equilibrium on the 
left is axisymmetric, and can 
exist in a star (giving a dipo- 
lar field on the surface) as 
well as inside a non-gravi- 
tating bubble. On the right, 
a more complex non-axisym- 
metric equilibrium, which can 
only exist in a star. [Blue 
shading represents surface of 
star.] 

by the (known) characteristic scale-length of the equilibrium. Unfortunately, we have few 
constraints on how much helicity a protostellar dynamo is likely to produce (or destroy). 

As mentioned above, there are a few mechanisms which cause a fossil field to evolve 
and which could also bring a core-dynamo field upwards to the surface. First, finite 
conductivity, which 'spreads out' the magnetic energy spatially; the timescale is ~ 10 10 yrs 
in a MS star. Second, buoyancy: imagine a magnetised region surrounded by a less or 
non-magnetised medium - since the field provides pressure without mass, the magnetised 
region must have a lower temperature in order to have the same total pressure and density 
as its surroundings. Heat diffuses inwards, buoyany rise on a timescale Tbuoy ~ P tk-h 
where j3 is the ratio of thermal to magnetic pressure (realistically > 10 6 ) and tk-h 
is the thermal timescale of the star so that Tbuoy 10 12 yrs. In neutron stars there is 
an equivalent mechanism where the role of temperature is replaced by electron fraction 
Y e . Finally, meridional circulation: a star in solid-body rotation cannot satisfy both the 
momentum equation and the heat equation everywhere (the von Zeipel paradox), and 
some additional circulation is required (von Zeipel 1924, Zahn 1992, Decressin et al. 
2009 and refs. therein.) This circulation takes place on the so-called Eddington-Sweet 
timescale te-s ~ Tk-h(-P /-^breakup) 2 where tk-h is the stellar thermal timescale. This 
is just possibly relevant for the fastest rotating Ap stars. 

In summary, we expect diffusive evolution to be negligible on the MS. In contrast, in 
neutron stars we expect the diffusive timescales to be rather shorter, in agreement with 
observational evidence of field decay. Some questions remain: why do not all A stars have 
these strong fields, and why is there such a large range in field strengths amongst the Ap 
stars - at what stage during formation is the Ap destiny of the star determined? 

3. Magnetism and Rotation in the Herbig Ae/Be stars 

Fossil field theory (see §2) predicts that fields are present throughout the MS and also 
earlier, in the pre-main sequence (pre-MS) phase. To explain the slow rotation of Ap/Bp 
stars, magnetic braking during the pre-MS has been proposed. 

To test both hypotheses, E. Alecian and collaborators performed a high-resolution 
spectropolarimetric survey of over 100 of the pre-MS progenitors of the A and B stars, the 
Herbig Ae/Be stars, discovering 7 magnetic stars among 128 observed (HD 200775, HD 
72106, V380 Ori, HD 190073, NGC 6611 601, NGC 2244 201, NGC 2264 83), implying 
that 5% of the Herbig Ae/Be stars are magnetic, as predicted by the fossil theory. 
Furthermore they performed a monitoring of the magnetic stars and were able to fit the 
temporal variations of their Stokes V profiles using the oblique rotator model (see an 
example for V380 Ori in Fig. They conclude that these stars host large-scale organised 
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magnetic fields with dipole strengths between 300 G and 3 kG, as predicted by the fossil 
theory (e.g. Alecian et al. 2008, Alecian et al. 2009, Wade et al., in prep.) 
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Figure 2. LSD V profiles of V380 Ori superimposed by the best oblique rotator model (dashed 
red line). The numbers close to the profiles are the rotation phase, and the small bars on the 
right of the profiles are the mean error bars in V. The profiles are sorted by increasing rotational 
phase, and the date of observation is indicated next to each profile (Alecian et al. 2009). 

Then they measured the rotation velocities projected on the line of sight (vsini) of all 
the stars of their sample, and they first compared the vsini distributions of the magnetic 
to the non-magnetic stars (Fig. [3]) . They find that the magnetic stars have been braked 
more than the non-magnetic ones, and that the braking must occur very early during 
the pre-MS evolution. Finally, they compared the v sin i distribution of the normal (i.e. 
non-magnetic and non-binary) HAeBe stars projected onto the zero-age main-sequence 
(ZAMS) to the normal A/B stars on the MS (Fig.lJ). They find that the two distributions 
are very similar and conclude that the normal HAeBe stars are expected to evolve towards 
the ZAMS with constant angular momentum (Alecian et al., in prep.) 



Figure 3. iisini his- 
tograms of the magnetic 
(left) and the non-mag- 
netic (right) HAeBe stars. 
The y-axes are graduated 
in percentage of stars on 
the left, and in number of 
stars on the right (Alecian 
et al., in prep.). 



Figure 4. iisini his- 
tograms of the normal 
HAeBe stars projected on 
the ZAMS (left) and the 
normal A/B stars (right). 
The y-axes are graduated 
in percentage of stars on 
the left, and in number of 
stars on the right (Alecian 
et al., in prep.). 
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4. Magnetic fields of O and B stars: measurements, statistics and 
evolution 

Recent measurements of the magnetic fields of early-type stars have been collected in 
a new catalogue of the magnetic fields of OB A stars ( Bychkov et al. 2009 1 . Based on 
the data from the catalogue together with the newest data, the statistical properties of 
an ensemble of the magnetic fields of OB stars were investigated. The rms longitudinal 
magnetic field B was used as a statistical measure of the mean field. This statistical 
experiment showed that this value depends weakly on the distribution of the random 
times of observations and on the structure of the field. The mean magnetic field averaged 
over the spectral subclasses has an unexpected jump between O and B spectral classes 
as can be seen in Fig. [S] (left panel) . 




[(B)=0.33.(B/1kG)" ,B2 



Spectral Class B, kG 

Figure 5. Left Panel: Mean magnetic fields of OB stars for different spectral subclasses. 
Right Panel: Magnetic field function for B stars. 

We have calculated the differential magnetic field function (MFF) F(B) for B stars, 
which is defined such that F(B)dB is the probability that a rms longitudinal mean field 
B lies in the interval (£, B + dB) (right panel of Fig. [Q). It is found that the MFF for 
B > 400 G can be approximated by the power function F(B) = A(B /lfcG) -7 with 
parameters A — 0.33 and 7 = 1.82. The MFF steeply decreases for B < 400 G - a result 



compatible with that of Auriere et al. (2007) who suggest that weak fields are destroyed 
by instabilities. It is also found (Kholtygin et al., in prep.) that the rms mean magnetic 
field B of the star can decrease by a factor 4-5 times during its evolution from ZAMS to 
TAMS in accordance with ?. 



5. The effect of atomic diffusion in stellar evolution 

Georges Michaud and his collaborators (O. Richard, J. Richer and M. Vick) are using 
abundance anomalies observed in AmFm and Horizontal Branch (HB) stars to constrain 
the hydrodynamics of stellar interiors. Their basic code includes all standard equations of 
stellar evolution to which they add 56 coupled differential equations to introduce atomic 
diffusion of all species included in the OPAL Rosseland opacities currently used to model 
stellar interiors in addition to Lithium, Berillium and Boron, which were calculated by the 
Montreal group. They use the original spectra that were used by OPAL to calculate their 
opacities. The dominant term in the transport equations involves the difference between 
gravitational and radiative accelerations. It is constantly recalculated for every mass shell 
during the evolution along with the Rosseland opacity so that the calculations take all 



composition changes consistently into account. For more detail see Richer et al. (1998) 



Turcotte et al.(1998) and Richard et al. (2001) 



Results indicate that atomic diffusion driven by radiative accelerations can explain 
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the abundance anomalies observed in AmFm and hot HB stars. In the absence of any 
competing hydrodynamical process, atomic diffusion actually leads to larger anomalies 
than observed. Two competing processes have now been investigated: turbulence mixing 
a given external mass, and mass-loss at a constant rate. Observations of Sirius A and 
of the Hyades star 68 Tau were found to be explained equally well by either of the two 
models. In the case of Sirius, 12 out of the 16 elements have their surface abundances well 
reproduced either by a mixed outer mass of about 10~ 5 solar mass or by a mass-loss rate 
of 10~ 13 solar mass per year. In the case of 68 Tau, 13 of the 14 elements are obtained 
within error bars. The mass-loss rate implies that after 10 8 years the matter at the surface 
largely comes from about 10 -5 solar mass below the surface which is also the mixed 
mass in the model with turbulence. That is where the important competition between 
gravitational and radiative accelerations occurs in both models. Observation of pre-MS 
abundance anomalies would favor mass loss as a competing mechanism. For a detailed 
description see |Richer et al. (2000)| for the turbulence model and Vick et al. (2009, in 
prep.) for the mass loss model. 

It is emphasized that there are large differences in the abundances determined by 
different observers and that a critical analysis of the abundance results is urgently needed 
to improve the situation and permit more rigorous constraints on the models. 




Figure 6. Concentration variations for all calculated species in two 2.50 solar mass models 
at 500 Myr; Left Panel with mass loss (10 _13 Mo/yr) as the competing process and Right 
Panel with turbulence as the competing process. An overabundance by a factor of 1.7 or more 
appears black in both cases while an underabundance by a factor of 0.5 or less appears white. 
The radial scale is linear in r. Horizontal lines indicate the mass of the spherical shell outside a 
certain radius (Am) labeled by log(Am/M). As is easily seen by comparing the two panels, the 
interior abundances are quite different in the outer 20 % of the radius of the star but, yet, both 
models have the same surface abundances which also agree reasonnably well with the anomalies 
observed on some AmFm stars. See Vick et al. (2009, in prep.) for details and a similar hgure 
at a different age. The anomalies caused by atomic diffusion affect the outer 30 % of the radius 
while nuclear reactions modify abundances of lighter elements within the inner 50 % of the radius 
leaving only a small buffer where the original abundances of the star have not been modified. 

However, even if the surface abundances obtained taking either mass loss or turbulence 
into account are very similar, the interior composition is very different - see Fig. [51 
The difference comes from the different nature of the flux in the two models. In the 
model with mass loss, the flux is advective and for each species approximately given by 
pc(r,t)(vD + vw) where p is the plasma density, c(r,t) is the local concentration of a 
given species at time t, vd is the diffusion velocity and vw the wind velocity. When the 
radiative acceleration is much larger than gravity, the vd is positive and adds up to vw 
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so that c(r, t) must be smaller for a given mass loss rate. In those regions of the envelope 
where radiative acceleration is largest on a given species, the concentration of that species 
is smallest in a mass loss solution. Heavy metals then tend to be underabundant between 
the surface convection zone and 10~ 5 solar mass below the surface whereas they are 
overabundant there in the model with turbulence. Helium also turns out to be distributed 
differently in the two models. This opens the interesting possibility of distinguishing 
between the two models by asterosismology. The outer 1O _5 M0, although a small mass, 
occupies approximately 20 % of the radius of the star and this is where the driving 
often occurs most efficiently for seismic waves. The two models should have significantly 
different asterosismic signatures. These have not yet been calculated in detail. 

On the HB, overabundances of Fe by a factor of order 50 are seen in all stars hotter than 
11,000K in the cluster M15 except for one star which has a relatively large v sini. This can 
be understood as the consequence of atomic diffusion driven by radiative accelerations as 
discussed in |Michaud et al. (2008)[ The appearance of abundance anomalies in HB stars, 
AmFm stars and HgMn stars is explained by their slow rotation if one uses a simple 
meridional circulation model which involves no adjustable parameters, as is discussed in 



Quievy et al.(2009) Since this model uses advective flows, it is tempting to assume that 



the concentration variations seen in the advective flow of the mass loss model would also 
apply to the meridional flow and this opens the possibility of a test by asterosismology 
of the nature of the process competing with atomic diffusion as one considers stars with 
faster rotation in which surface anomalies progressively disappear. 

Diffusion in magnetic stellar atmospheres Most magnetic CP stars show abun- 
dance variations that are the signature of separation in the atmosphere itself. The mod- 
elling of transport processes in the atmosphere is much more demanding than in stellar 
interiors since the magnetic fields necessitate 3-D modelling for a complete description, 
which in principle requires 3-D radiative transfer and particle transport. It should lead 
to the appearance of abundance spots within the observed magnetic structures. 

Results of time dependent diffusion of strontium in magnetic atmospheres were pre- 
sented by G. Alecian. In collaboration with M. Stift he solved in detail the particle 
transport separately in the presence of horizontal and vertical magnetic fields. The two 
solutions are very different. The accumulation can be very short, of the order of days 
in the high atmosphere. They suggest that as a first approximation, one may consider 
transport in the vertical direction only, since the distances are shorter; the spots are then 
created by the effect of the magnetic field on the vertical diffusion velocity. These results 
are very encouraging but many steps are still required to obtain a map of abundance 
anomalies at the surface of magnetic CP stars. 

On the other hand, it was suggested by J. Portnoy and R. Steinitz using qualitative 
arguments that magnetic field gradients could lead to an additional term in the particle 
transport equation, perhaps explaining the over-abundance of heavier elements, and that 
dielectronic recombinations could be potentially important for the ionization equilibria 
in magnetic atmospheres. 



6. Towards a coherent picture of internal transport in CP stars 

The study of asteroseismology and powerful ground-based instrumentation dedicated 
to stellar physics is developing strongly (CoRoT, KEPLER, ESPaDOnS, etc.) generating 
tight constraints on the internal structure, dynamics, and magnetism of stars. For this 
reason stellar models are needed that include transport processes from the birth of stars 
to their death. A coherent picture of the dynamics of stellar radiative zones, where non- 
standard chemical mixing occurs, is thus required (cf. Zahn 1992, Meynet & Maeder 
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2000). A complex transport, the rotational transport, which involves several mechanisms, 
takes place in these regions. 

First, rotation, structural adjustments and angular momentum losses at the surface 
induce 'baroclinic' or 'meridional circulation' (see $2j Maeder & Zahn 1998, Mathis & 
Zahn 2004, 2005) which acts simultaneously to transport angular momentum, chemicals 
and magnetic field by advection as well as inducing differential rotation. Next, differen- 
tial rotation induces hydrodynamical turbulence via shear, baroclinic, multidiffusive and 
other instabilities, which of course in turn damps the differential rotation - as happens 
in the terrestrial atmosphere (Zahn 1983; Mathis et al. 2004 and refs. therein). Then, 
(differential) rotation interacts with fossil magnetic fields. In this case, the mean secular 
torque of the Lorentz force (Gough & Mclntyre 1998, Mathis & Zahn 2005, Brun & Zahn 
2006) and magnetohydrodynamical instabilities such as Tayler and multidiffusive mag- 
netic instabilities (Tayler 1973, Spruit 1999, Menou et al. 2004) can modify the transport 
of angular momentum and chemicals. Eventually, a dynamo in stably-stratified stellar 
radiation zone is possible due to the non-linear interaction of instability-velocity field 
and instability-magnetic field but its detection in numerical simulations is still under 
exploration (Spruit 2002, Braithwaite 2006, Zahn et al. 2007). 

Finally, internal gravity waves are excited at the boundaries with convective zones. 
They propagate through radiative regions where they extract or deposit angular momen- 
tum at the location where they are damped, leading to a modification of the angular 
velocity profile and consequently of the chemical distribution (Goldreich & Nicholson 
1989; Talon & Charbonnel 2005). Because of their frequencies, they can be strongly 
modified by both the Coriolis acceleration (Mathis et al. 2008, Mathis 2009) and the 
Lorentz force (Kumar et al. 1999). 

If the star is one component of a close binary system, tidal processes modify the angular 
momentum transport via the external torque of the equilibrium tide and the dynamical 
tide, i.e. the internal gravity waves which are excited by the tidal potential (Zahn 1977). 
Like external winds, tides can strongly modify the internal angular velocity distribution 
and thus the magnetic behaviour and mixing. 

The rotational transport and various interactions are summarized in Fig. [7] These 
mechanisms now have to be applied to massive peculiar stars, to study their rotational 
history and their potential equilibrium state such as the magnetohydrostatic equilibrium 
(see Moss 1977, Duez et al. 2009, and §2). 



Figure 7. The com- 
plete rotational trans- 
port in a stellar radia- 
tion zone. 
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